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ABSTRACT 

Many of the compact extragalactic radio sources that are used as fiducial points to define 
the celestial reference frame are known to have proper motions detectable with long-term geode- 
tic/ astrometric Very Long Baseline Interferometry (VLBI) measurements. These changes can 
be as high as several hundred micro-arcseconds per year for certain objects. When imaged with 
VLBI at milli-arcsecond (mas) angular resolution, these sources (radio-loud active galactic nu- 
clei) typically show structures dominated by a compact, often unresolved "core" and a one-sided 
"jet". 

The positional instability of compact radio sources is believed to be connected with changes in 
their brightness distribution structure. For the first time, we test this assumption in a statistical 
sense on a large sample, not only for individual objects. We investigate a sample of 62 radio 
sources for which reliable long-term time series of astrometric positions as well as detailed 8-GHz 
VLBI brightness distribution models are available. We compare the characteristic direction of 
their extended jet structure and the direction of their apparent proper motion. 

We present our data and analysis method, and conclude that there is indeed a correlation 
between the two characteristic directions. However, there are cases where the ~1 — 10-mas 
scale VLBI jet directions are significantly misaligned with respect to the apparent proper motion 
direction. 

Subject headings: radio continuum: galaxies — galaxies: active — quasars: general — techniques: 
interferometric — reference systems — BL Lacertae objects: individual (OJ 287) 
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1. Introduction 

The International Celestial Reference Frame 
(ICRF) is a quasi-inertial reference frame defined 
through the accurate radio positions of extragalac- 
tic sources (active nuclei of quasars, radio galaxies, 
and BL Lac: objects) distributed over the sky. Its 
first realization (ICRFl: lMa et al.ll998l ) was based 
on the positions of 212 defining radio sources, 
using Very Long Baseline Interferometry (VLBI) 
measurements made between 1979 and 1995. 

The continuous accumulation of new high- 
quality VLBI data, and technical, observing and 
data analysis impro vements made i t possible to 
define the ICRF2 (|Fev et all l2009h which was 



significant apparent proper motion on a scale 
up to 100 micro-arc seconds (/xas) per year (e.g. 
Feissel- Vernier! 120031 ). These systematic linear 



adopted by the International Astronomical Union 
(IAU) in 2009 as the current official realization of 
the celestial reference system. There are 295 defin- 
ing sources selected for ICRF2, based on their po- 
sitional stability and compact radio structure, as 
deduced from nearly three decades of VLBI mon- 
itoring. The ICRF2 catalog lists the coordinates 
of a total of 3414 extragalactic radio sources. The 
new solution is aligned with the ICRF1 to ensure 
continuity. The new reference frame contains the 
positions of more than 5 times as many sources, 
better distributed over the northern and south- 
ern hemispheres. Its accuracy is at least 5 times 
better than that of the ICRF1. 

Due to the frequent observations in the frame- 
work of geodetic/astrometric VLBI programs, 
the positional accuracy of many of the reference 
sources is better than ~0.2 milli-arcseconds (mas) 
in both equatorial coordinates. It has long been 
known, however, that several of these sources show 
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or non-linear variations in the position of some 
reference sources can lead to inaccuracies in the 
estimates of other p arameters (e.g. Earth ori- 
entation parameters: iMacMillan fc Mai 120071 and 
references therein). 

The typical angular scale of ~1 — 10 mas of the 
emitting regions in these radio sources that can be 
studied with VLBI imaging at centimeter wave- 
lengths is significantly larger than the formal pre- 
cision of the position determination. Variations 
in the brightness distribution of quasars and radio 
galaxies are not rare events. It is a general as- 
sumption that in most cases the apparent proper 
motions are mainly due to this intrinsic structural 
evolution, e.g. motions of relativistic jet com- 
ponents away from the innermost regions of ac- 
tive galactic nuclei (AGNs). This hypothesis was 
confirmed for some individual objects. For exam- 
ple, a detailed analysis of 0923+392 (4C +39.25) 
showed that the apparent proper moti on of the 
sourc e is related to structure variations ( Fev et all 
19971 ). The apparent proper motion of 2200+420 
(BL Lac) was found to significantly decrease after 
structural corrections derived from source m aps 
were applied to the data set (jCharlod 12002). In 
some cases, the proper motions of the q uasars 
are n ot easily modeled by a linear function (jTitovl 



2007). These investigations implied that the posi- 



tional instabilities can be explained by examining 
the variations of the source structure in these indi- 
vidual cases. However, the general hypothesis has 
not been evaluated on a larger sample yet. 

If this model can be applied for most of the 
sources, then we can expect some correlation be- 
tween the direction of the position change, and 
the typical direction of the source "core-jet" struc- 
ture in the sky. The mas-scale radio structure of 
AGNs is dominantly one-sided. This is a conse- 
quence of relativistic beaming which enhances the 
radiation originating from the relativistic plasma 
outflow pointing towards the observer, close to 
the line of sight. At the same time, the reced- 
ing side of the intrinsically symmetric jet is de- 
boosted and becomes practically invisible. The 
compact, often unresolved, flat-spectrum "core" 
one can usually see in VLBI images is in fact 
the base of the inner jet, characterized by an 
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optical depth of r = 1. The absolute position 
of this "core" therefo re depends on the observ- 
ing frequency (e.g. iBlandford fc Konigll 19791 
Lobanovl Il998l iKovalev et all 120081) . Even at a 
fixed frequency, positional changes are expected 
with time due to changes in the intrinsic physical 
(e.g. jet Lorentz factor, electron density) or ge- 
ometric (viewing angle, jet precession) conditions 
within the innermost few parsecs of the AGN en- 
vironments. These changes manifest themselves in 
brightness, spectral or polarisation variations, ap- 
pearance, disappearance, enhancement or fading 
of optically thin VLBI "jet" components. These 
phenomena are detectable with high-resolution 
VLBI imaging typically over time scales of years 
or longer. Propagation effects (e.g. interstel- 
lar scattering) extrinsic to the compact sources 
could also lead to apparent structural or bright- 
ness modifications, down to s ub-daily time scales 
(e.g. IWagner fc Witzellll995l) . The changes could 
in principle affect the absolute position of the ra- 
dio centroid of a radio-loud AGN which serves as a 
reference point in astrometric and geodetic VLBI 
observations. 

Some systematic effects could contribute to the 
total magnitude of the apparent proper motion 
as estimated from geodetic VLBI data. One of 
them is the secular aberration drift caused by 
the Galactocentric acceleration of the Solar sys- 
tem barycenter. This acceleration a — V 2 /R 
(where V is the orbital velocity due to rotation 
of the Solar system around the Galactic center, 
and R is the Galactocentric distance) estimated as 
2.7 x !CT 13 kms- 2 (V = 254km s _ \ R = 8.4kpc, 



Reid et al. 2009) results in a secular aberration 



drift equal to 5.4 /las/year. The effect of noncir- 
cu lar motion was found neglig ible (<1 ^as/year) 
by lKopeikin fc Makarovl (|2006h . 

iTitov et"aL ( 2011 ) detected the secular aberra- 
tion drift in radio source proper motions as derived 
from 30 years of the geodetic VLBI observations. 
They found 6.4±1.5 /zas/year, which is consistent 
with the theoretical prediction. Thus, this effect 
should be removed from the apparent proper mo- 
tion of individual radio sources for a fair compar- 
ison with astrophysical data. 

In the present work, we compared the charac- 
teristic position angle of the core-jet structure of 
62 radio sources derived from their multi-epoch 
VLBI structure maps with the direction of their 



apparent proper motion based on the coordinate 
time series of regular geodetic VLBI observations. 
The individual proper motions have been cor- 
rected for the secular aberration due to the ro- 
tation of the Solar system barycenter around the 
center of the Galaxy. The data sets and the data 
analysis are described in § [2l The results are pre- 
sented in § [3l Implications of our results are dis- 
cussed in § [4j 

2. Data 

2.1. Description of the apparent proper 
motion data 

To derive the apparent proper motion of radio 
sources, we analyzed more than 5,000 geodetic 
VLBI sessions of the Internation al VLBI Service 
for G eodesy and Astrometry (IVS. ISchluter fc Behrendl 
20071 ) data base since 1979 using the Calc 10.0/Solve 2010.05.21 
geodetic VLBI analysis software package, devel- 
oped and maintained at the NASA Goddard Space 
Flight Center. A more detailed de scription of 
the an alysis configuration is given in ITitov et al.l 
([20 111 ). Earth orientation parameters, station and 
radio source coordinates were estimated once per 
session. Tropospheric zenith delays and gradients 
were estimated at a higher rate. We used a loose 
constraint t o tie the celest ial reference frame to 
the ICRF2 (|Fev et al.ll2009h in order to prevent a 
spurious net rotation of a set of radio sources that 
could result in a fake orientation of the proper 
motion. 

The above analysis configuration allowed us to 
obtain coordinate time series for all the observed 
radio sources. In these series, data points result- 
ing from less than three good observations within a 
session were removed and outliers were eliminated 
so that the reduced \ 2 is reasonably close to unity. 
Proper motions were then estimated for time se- 
ries containing at least 10 points and longer than 
10 years by weighted least-squares with weights 
taken as the inverse of the squared formal error. 
The final sample contains the proper motions of 
593 sources. 

The proper motions were freed from the secu- 
lar aberration drift caused by the rotation of the 
Sun around the Galactic cent er. The contribu- 
tion is given by Eqs. (3)-(4) of ITitov et ail (|201ll ) 
wherein the coefficients di correspond to a vec- 
tor of amplitude of 6.4 /xas/year directed toward 
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Description of the VLBI structure 
data 



In order to determine the characteristic posi- 
tion angle of the core-jet structures in a large 
number of radio sources, and to monitor the vari- 
ability of the typical je t directions , we us ed the 
database compiled by IPiner et~aL ( 2007 . here- 
after P07). This data set is based on the US 
Naval Observatory's Radio Reference Frame Im- 
age Database, an ongoing program t o image ra- 
dio re ference frame sources regularly. IPiner et al 



(l2007h used the data series of 77 sources obtained 
at 8 GHz. They fitted Gaussian brightness distri- 
bution model components interactively to the vis- 
ibility data on a source-b y-source basis with th e 
Caltech Difmap software ([Shepherd et al. 1 11994 . 
Their final database (see Table 3 in P07) includes 
the main parameters of the fitted components (e.g. 
the polar coordinates of the centers, the sizes and 
orientations of the major and minor axes, and the 
flux densities of the Gaussians) for each source 
at different epochs. All of the 77 sources in P07 
can be found in the proper motion database as 
well. We note that 47 of them are ICRF2 defining 
sources, and thus are among the most stable and 
less prominently extended objects. 

2.3. Data analysis 

For the selected sources, the P07 database pro- 
vides model brightness distributions at 3-19 dif- 
ferent epochs over the years 1994-1998. We fitted 
a straight line to the positions of the components 
that describe a specific source using a weighted 
least-squares method. A position angle was cal- 
culated independently at each epoch. The fit was 
weighted by the flux densities of the individual 
emission features. Model components with ID=99 
(Table 3 in P07) were not considered in our anal- 
ysis, since these components could not be directly 
identified with model components that appear at 
other epochs in P07. The final characteristic core- 
jet direction (cj>) and its uncertainty (er^) were 
computed as the average and the standard devi- 
ation of the position angles obtained at different 
epochs. Fig. [1] illustrates our fitting method for 
0202+149. 



The apparent proper motion values in right as- 
cension and declination, and their uncertainties 
were transformed into polar coordinates, yielding 
four parameters for each source: the radial proper 
motion (r) and the position angle (tp) with their 
uncertainties (ay, cy,). Sources with negligible ap- 
parent proper motions (r/oy < 1) were discarded 
from our sample and were not used in the further 
analysis. The derived parameters of the final se- 
lection of 62 sources are summarized in Table [T] 

3. Results 

The variations present in the apparent posi- 
tions of AGNs are believed to be mainly caused 
by changes in the brightness distribution of the 
sources. According to this general picture, the di- 
rection of the apparent proper motion of a specific 
source should show correlation with the typical di- 
rection of the source core-jet structure in the sky. 
The multi-epoch VLBI structure maps offer a good 
opportunity to check whether the typical core-jet 
directions are stable in time for our sample. Both 
the real variations in the core-jet direction and 
the uncertainties of the fitted Gaussian models can 
contribute to the uncertainties of the derived posi- 
tion angles. For most of our sources, these values 
are relatively low (<20°), implying that the role 
of real variations may be negligible. Therefore in 
the further analysis we assumed that the derived 
directions are valid on longer time scales as well. 
The time span of the geodetic VLBI observations 
is up to five times longer than that of the structure 
maps analyzed in P07. 

To verify our assumption, we studied an ex- 
treme example, the source OJ287 (0851+202) 
in more detail. The blazar OJ 287 is known to 
have a rapidly changing inner rad io jet structure 
(e.g. iTatevama fc Kingham 2004 . and references 
therein) of which the kinematics could be revealed 
by frequent VLBI monitoring. We investigated 
the mas-scale jet direction in OJ 287, as defined 
by the component positions in 8-GHz VLBI imag- 
ing data, over a time interval much longer than 
covered by P07. To this end, we collected data 
from the literature (1985-1988, 1990-1996), and 
analyzed archival geodetic/astrometric VLBI vis- 
ibility data (1994-2008). The total time range 
spanned is 24 years, comparable to the time range 
of the data from which the proper motion is esti- 
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mated (§13}. 

A total of 59 epochs of 8-GHz Very Long Base- 
line Array (VLBA) data sets were taken from the 
VLBI database of compact radio sourcefl from 
1994 to 2008. In comparison, P07 used data from 
15 of these epochs, between 1994 and 1998. With 
the Difmap software, we model-fitted the cali- 
brated interferometric visibilities at 40 additional 
epochs where sufficient data were available, adopt- 
ing the same method used by P07. In addition, to 
extend the covered t i me in terval backwards, we 
refer to lVicente et al. ( 19961 ) who analyzed 8- GHz 
geodetic VLBI observations at 10 epochs (jet com- 
ponent was detected in nine epochs) between 1985 
and 1988. They explain the observed radio struc- 
ture with a helical jet arising from a supermas- 
siv e binary black hole sy stem. A similar study 
bv lTatevama et ah ( 19991) used 27 epochs between 
1990 and 1996. 

Fig. H shows the characteristic jet position an- 
gle of OJ 287 as a function of time. There is in- 
deed a clear systematic change in this angle, quite 
consistent with the predictions of the ballistic pre - 
cessing jet model of iTatevama fc Kinghaml ( 2004 ). 
However, from the point of view of the present 
study, it is important to note that the jet always 
extends nearly to the west-southwest (position an- 
gles between about -80° and -135°). Within 
the uncertainties, the characteristic core-jet direc- 
tion (<p) for OJ 287 we determined from the P07 
data (§ 12. 3p agrees with the average value over 24 
years. In other words, our method of using only 
the 5-year P07 data is justified since it does not in- 
troduce a significant bias in the characteristic jet 
direction, even in the case of this peculiar radio 
source. 

The changes in the brightness distribution of 
the radio sources could cause apparent positional 
variations in the opposite directions along the 
core-jet axis. In order to define a modulo 90° 
differential angle between the the two directions, 
we computed a <p angle for each object using the 
following formulae: (i) <p = <p, if the difference 
between 4> and V is <90°; (u) 4 = 4> + 180°, if 
the angle between cp and ip is >90° and (f> < 0°; 
(Hi) <p = (j) — 180°, if the angle between <p and 
ip is >90° and (p > 0°. The differential angle 
(A$ = \ip — <f) \) is therefore the smallest angle 



between the direction of the apparent proper mo- 
tion and the axis defined by the core-jet structure. 
The derived A$ parameters and their cta* uncer- 
tainties are listed in the last column of Table [TJ 

3.1. General association between the char- 
acteristic directions 

Since ip and <p are circular variables, the anal- 
ysis of the possible association between them re- 
quires special methods. In order to test the hy- 
pothesis whether the two variables are indepen- 
dent, we performed a g eneral t est, fo llowing the 
procedures described in iFisher (1995). As a first 
step of this analyis, ip and <p were re-ordered 
to get data pairs of (ip*,^)*),...,(ip*,(f)*), where 
ip*,...,ipn are in cyclic order, and Si,...,S n were 
determined as the rank of the corresponding (p * 
values. Then we calculated 

n n 
i=Q j=Q 

where 

T{j = n ■ min(i, Sj) — i - Sj, i,j = 1, n. 

The total number of radio sources in the analy- 
sis was n = 62. The null hypothesis that ip and 
<p> are independent can be rejected if 7^ is too 
large. In order to determine a critical value for 
7^, we perfo rmed a randomisation test (see also 



'http: //astrogeo . org/vlbi_images/ 



Fished [19951) : (a) we added uniformly-distributed 
random numbers in the range of 0° and 360° to 
the position angles derived for the core-jet struc- 
tures (4> ) and then re-calculated the value of 7^; 
(b) this treatment was repeated 1000 times, and 
the 7^ values obtained were compared with the 
original one. We found that the value of 7^ for 
our original sample is quite high: the test showed 
that the null hypothesis, that the two variables are 
independent, can be rejected at a level of 98.3%. 

The distribution of the A$ values is presented 
in Fig. [3] (left). The histogram shows a peak at 
low differential angles, which could be expected 
if there is correlation between the two directions. 
Using the results of the previous randomisation 
test, we also compiled the simulated distribution 
of A$ (Fig. El left) and A$/cr A $ (Fig. S right) 
by averaging the results of the 1000 different runs 
in each bin. As Fig. [3] indicates, both distribu- 
tions show an excess compared to the simulated 
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distributions at their lower end. This finding also 
implies that there is an association between the 
two characteristic directions. For some sources 
this association seems to be very strong. The best 
examples were selected based on the following cri- 
teria: the differential angle between the two typ- 
ical directions is A<£> < 5°, and A$/cta$ < 1. 
Two objects fall into this category, 0919-260 and 
1156+295. 

On the other hand, there are a couple of sources 
where the angle between the direction of the ap- 
parent proper motion and the axis defined by 
the core-jet structure is quite large. 0851+202 
(OJ 287) is especially interesting from this point 
of view, since this object forms a clear outlier at 
A$/(ja* = 6.4 in Fig. [3] (right) showing that its 
nearly orthogonal misalignment (see Table [1]) is 
highly significant. 

We note that the omission of the correction for 
the secular aberration in the calculation of proper 
motion data would not alter our main conclusion 
concerning the association between the two char- 
acteristic directions. To test this, we repeated the 
analysis using a data set derived without this cor- 
rection. In that case, the null hypothesis could 
be rejected at a level of 99.5%. Before 1990, the 
general deficiency of the VLBI networks, including 
the number of observed sources and observing an- 
tennas per sessi on, makes the VLBI products less 
reliable (see. e.g.lGontier et al.ll200lHMalkinll2004t 



Feissel- Vernier et al.l 20041 lLambert fc Gontier 
2009, who reported interesting statistical results 
and remarks about the VLBI evolution over the 
last two decades). In order to check the possible 
effect of these deficiencies on our findings, we re- 
peated our analysis using a recompilation of the 
proper motion database without the geodetic data 
measured before 1990. We obtained a significance 
level of 96.6% for the rejection of the null hypothe- 
sis. It is in good agreement with the value derived 
utilizing the original data set. 

4. Discussion 

Based on our sample of 62 radio-loud AGNs, 
there is correlation between the direction of their 
apparent proper motion and the typical direction 
of their core-jet structure seen in cm-wavelength 
VLBI images. However, this does not necessarily 
mean that the proper motion is only influenced by 



changes in the source brightness structure. There 
exist individual objects where the apparent proper 
motion direction and the characteristic core-jet 
structure direction differ significantly. 

In fact different angular scales are probed by 
the 8-GHz VLBI imaging and the astrometric po- 
sition measurements. While the typical angular 
size of the core-jet structures we investigated is 
~1 — 10 mas, the measured proper motion values 
are at most in the ~0.01 — 0.1 mas/year range. 
The characteristic inner jet direction could be dif- 
ferent within the "core" that remains unresolved 
in the VLBI images at 8 GHz. In the exten- 
sive 15-GHz Monitoring Of Jets in Active Galac- 
tic Nuclei (MOJAV E) survey with the VLBA, 
Homan et al.l ( 20091) found that jet components of- 
ten change direction after ejection. The motion 
of about half of the 203 components they studied 
is significantly non-radial. These misalignments 
are typic ally within 30° but can reach up to 80° 
(Fig. 3 of iHoman et~aLll2009h . Misalignments be- 
tween the apparent jet position angles in com- 
pact extragalactic radio sources are well known 
at larger angular scales as well. The distribu- 
tion of the misalignments of the mas-scale struc- 
tures probed by VLBI, and the arcsecond-scale 
jets seen in the images made by the connected- 
element Very Large Array (VLA) interferometer 
sho ws a bimodal form, with peaks close to 0° a nd 
90° (|Conwav fc Murphvl [19931 lAppl et al.lll996l ). 

In our study, the nearly orthogonal misalign- 
ment found in 0851+202 (OJ287) could pos- 
sibly be explained with a similar effect, on a 
smaller angular scale. OJ 287 is a source known 
to ha ve a precessing jet (jTatevama &: Kinghaml 
2004, and references therein). Interestingly, 



the highest-resol ution VLBI im age available for 
O J 287 (86 GHz, iLee et alJEooi l shows a hint on 
a ~100 /xas-scale weak extension to the compact 
unresolved core, close to the direction of the ap- 
parent proper motion measured for this source 
(-15.6°; Table [J). 

Proper motion directions could be more favor- 
ably compared with VLBI images made at an or- 
der of magnitude higher resolution. This can be 
achieved by decreasing the observing wavelength 

(A) and/or by using Space VLBI (SVLBI) to in- 
crease the length of the interferometer baselines 

(B) . (The angular resolution of an interferome- 
ter array is proportional to X/B.) For example, a 
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recent, and currently t he most extend ed, global 



86-GHz VLBI survey (|Lee et al.l 120081 contains 



single-epoch information on structural properties 
of more than 100 extragalactic radio sources at 
sub-mas scale. A practical complication is that 
the extended optically thin steep-spectrum fea- 
tures fade away at high frequencies. Therefore it 
may prove difficult to use the mm- VLBI images 
which are more dominated by a single "core" com- 
ponent to define a characteristic jet direction for a 
large sample of objects. The technique of SVLBI 
offers an increase in resolution in another way. 
It involves an antenna orbiting around the Earth 
which is co-observing with a ground-based radio 
telescope network. The recently canceled second- 
generation Ja panese SVLBI mission ASTRO-G 
(lTsuboil l2008) would have provided ground-space 



baselines that exceed the Earth diameter by a fac- 
tor of ^3. 

It is possible that systematic effects other 
than changes in the source brightness distribution 
structure also affect the measured apparent proper 
motions. The estimated dipole and quadrupole 
harmonics in the pattern of the vector field of the 
prop er motions in the sky (|MacMillan l2005HTitovl 



2009, and references therein ) could in principle 



be interpreted in several ways: as a consequence 
of the galactocentric rotation of the Solar system 
(note that our proper motion data have already 
been corrected for the estimated amount of this 
systematic effect); as a result of an anisotropic 
Hubble expansion; as an indication of the pri- 



verse (Kristian & Sachs 1966: 


3 vne et al. 


1996; 


Gwinn et al.lll997: Sovers et al.l 


1998; MacMillan 


2005; Titovl 20091). Observationally establishing 



the magnitude of the all-sky systematic effects is 
complicated by the fact that the reference source 
distribution is uneven, especially in the south- 
ern hemisphere. For instance, correlation coeffi- 
cients between the first and second degree spher- 
ical harmonics could reach 0.8-0.9 due to the 
paucity of the radio sources below —40° decli- 
nation (|Titov fc Malkinll2009h . 



Gontier et al.l (|200ll ) suggest that at low de- 



clinations (| S | < 20°), tropospheric mismodelling 
could lead to more uncertain position estimates, 
preferentially in the north-south direction. We 
could not detect any declination-dependent effect 
in the values of A$, maybe due to the small size 



of our sample. 

In our analysis presented above, we assumed 
that the apparent proper motion of the sources 
can be represented by a linear change as a func- 
tion of time. The VLBI position time series we 
used span nearly a quarter of a century for many 
of the objects. While the linear approximation 
is valid for most of the sources, the problem in 
general is clearly more complex. Physically, many 
of the sources could have undergone multiple jet 
component ejections over the time range stud- 
ied. Although the position angle of the jet tra- 
jectories (the path of the subsequent components) 
could remain similar for a given source, these pat- 
terns in the jet flow may change their apparent 
speed. Indeed, there are a few AGNs known to 



show non-linear pr oper motions (e.g. lTitovl l2007: 
Zharov et al.l 120091 ) . Coordinate time series for 
the most extreme cases are available in the ICRF2 
document (Chapter 4, Fev et al. 20091) . Therefore 
the linear fit applied here is not necessarily the 
best representation of their coordinate changes. 
This fact is reflected in the higher formal errors 
obtained for the fitted linear proper motion val- 
ues for these sources. However, here we analysed 
the directions, and used the fitted absolute val- 
ues of the apparent proper motion only to select 
the objects that show significant motion. Future 
extensions of our work could look for directional 
changes as well, since the apparent proper motions 
could follow curved paths. It would also be inter- 
esting to find a relation between the magnitude 
of the detected proper motions and the structural 
properties of AGNs. 

5. Conclusions 

We studied a sample of 62 compact extragalac- 
tic radio sources for which reliable long-term VLBI 
time series of astrometric positions as well as de- 
tailed 8- GHz VLB I brightness distribution models 
(|Piner et al.ll2007t) were available (Table [T]). We 
compared the direction of their apparent proper 
motion and the characteristic direction of their 
~1 — 10 mas-scale radio jet structure. Earlier, for 
a number of individual sources, the observed as- 
trometric proper motions have been found to be 
in a good agreement with the apparent motion of 
bright jet component s in the mas-scale radio struc - 
ture (e.g. 4C +39.25. iFev et al.lll997l : lTitovll2007t ). 
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However, the assumption that the apparent proper 
motions are closely related to changes in the radio 
brightness distribution has not yet been tested on 
a large sample of sources. With a statistical anal- 
ysis, we found that there is a general correlation 
between the direction of proper motion and the 
characteristic direction of core-jet sources. How- 
ever, there are some individual objects where this 
relationship does not hold. In particular, a nearly 
orthogonal misalignment was found in the case of 
O J 287, a blazar with a known precessing jet. Our 
results imply that the apparent proper motion is 
in general influenced by the changes in the source 
brightness structure. This could be better veri- 
fied in the future with higher-resolution mm-VLBI 
and/or Space VLBI imaging data. 

We are grateful for the referee's insightful and 
constructive comments which certainly improved 
the presentation and clarity of our work. This 
work was partly supported by the Hungarian Sci- 
entific Research Fund (OTKA K72515). The au- 
thors thank Lajos G. Balazs and Csaba Kiss for 
useful discussions. This research made use of the 
astrogeo.org database of brightness distributions, 
correlated flux densities and images of compact 
radio sources produced with VLBI, maintained by 
Leonid Petrov. This study also made use of the 
International VLBI Service for Geodesy and As- 
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Table 1 

The properties of the 62 selected radio sources. 
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Table 1 — Continued 



Basic properties Structure maps Geodetic VLBI observations 
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Note. — Col.(l): Name of the radio source derived from B1950 equatorial coordinates. Col. (2): Optical class of the source. 
B — BL Lac object, Q — quasar, G — galaxy, HP refers to high polarization. Th e class information was taken from 
iVeron-Cettv & Veronl i2003fl . Col. (3): Redshift values ( adopted from lPiner et al.|[2007h . Col. (4): Shows whether the specific 
object is a defining source of the ICRF1/ICRF2. D— 1 and 2 mark ICRF1 and ICRF2 defining sources, respectively. D— 3 
indicates sources that arc defining objects in both ICRF1 and ICRF2. Col. (5): Position angle of the final characteristic core- 
jet direction (0) and its uncertainty (c^). Position angle is measured in degrees from north through cast. Col. (6): Number 
of observations used to compute 0. Col. (7): Radial component of the apparent proper motion (r) and its uncertainty (oy) in 
/^as/year. Col. (8): Position angle of the apparent proper motion {ip) and its uncertainty (c^) in degrees. Col. (9): Number 
of observations used to compute ip. Col. (10): The smallest angle between the direction of the apparent proper motion and 
the characteristic core-jet direction in degrees. 
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CO 

03 
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Fig. 1. — Structure maps represented by fitted Gaussian models for 0202+149 at four different epochs, taken 
from P07. The different components of the radio source are denoted by circles, whose diameters correspond 
to the full width at half maximum (FWHM) size of the fitted Gaussians. The dashed lines show the results of 
our line fittings, while the grey areas correspond to the 0±cj0 region (derived from the position observations 
at 13 different epochs). 
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Fig. 2. — The characteristic position angle of 0851+202 (OJ287) core-jet structure as a function of time 
(gray diamonds) , as derived from 8-GHz VLBI imaging data taken between 1986 and 2008. The data from 
Piner et al. (2007) used for the analysis presented in this paper (filled circles) cover a much shorter time 



range (1994-1998). The source is known to have a precessing radio jet which explains the clear directional 
variations seen in this plot. Even for this source with a high ly and systematically variable structure, the 
average position angle determined from the iPiner et al . (2007) data (-99?5 ± 8?3) and from the data taken 
over more than two decades (— 103?6 ± 12?3) are similar within the uncertainties. 
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Fig. 3. — Distribution of the A$ (left panel) and the (right panel) values for the 62 radio sources 
analyzed (histograms). Average values obtained from 1000 similar but simulated random samples are shown 
with vertical bars whose size indicates the standard deviation. 
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